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We investigate the influence of a homogeneus and constant strong external magnetic field on 
the heavy-meson spectrum. Quarkonium states cc and bb are described within a non-relativistic 
framework and by means of a suitable potential model based on the Cornell parametrization. In 
particular, in this work we propose a model which takes into account the possible anisotropies 
emerging at the level of the static quark-antiquark potential, as observed in recent lattice studies. 
The investigation is perfomed both with and without taking into account the anisotropy of the static 
potential, in order to better clarify its effects. 
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I. INTRODUCTION 

In recent times there has been a great interest regard¬ 
ing the physics of strongly interacting matter in the pres¬ 
ence of strong external magnetic fields^i.e. such that^ 
eB ~ or larger (see, e.g.. Refs. [l|, 0 for recent re¬ 
views) . This topic could be relevant to the study of some 
dense astroph;gical objects, like magnetars Q, and for 
cosmology (j, 1^. However, the main interest was trig¬ 
gered by the fact that magnetic fields of this order of 
magnitude can be created in a laboratory, in particular 
in heavy-ion collisions , when two relativistic heavy 
ions collide with a non-zero impact parameter, producing 
a huge field in the collision region. For example, one can 
reach \e\B ~ 0.2 — 0.3 GeV^ in Pb-|-Pb collisions at the 
Large Hadron Collider (LHC). 

Such huge magnetic fields are produced in the very 
early stages of the collision. It is still not clear how long 
and to what extent they survive the thermalization pro¬ 
cess of the fireball created after the collision. Therefore, 
while various theoretical investigations, based both on 
model studies and on Lattice QCD simulations (LQCD), 
have predicted many interesting phenomena affecting the 
properties of strongly interacting matter in the presence 
of strong magnetic backgrounds, it is still uncertain to 
what extent such phenomena will be detectable in heavy 
ion experiments. 

In this perspective, effects regarding the physics of 
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^ For reference, a magnetic field of the order of Tesla corre¬ 
sponds to eB ~ 3.3 zz 0.06 GeV^. 


heavy flavors are of particular interest, since they are 
more sensitive to the conditions taking place in the early 
stages of the collision. Various studies have approached 
the issue of quarkonia spectra and production rates in 
the presence of magnetic backgrounds [iM3- Many in¬ 
teresting phenomena have been predicted, including the 
emergence of magnetic field induced mixings between dif¬ 
ferent states and of production anisotropies with respect 
to the collision plane. 

The starting point for most of these investigations is 
the coupling of the magnetic field with electric charges 
and magnetic moments carried by the valence quarks. 
However, the magnetic field is known to induce impor¬ 
tant modifications also at a non-perturbative level and in 
the gluon sector: a natural question is whether that can 
change the picture both quantitatively and qualitatively. 

A very interesting phenomenon in this respect is the 
magnetic field induced modification of the static quark- 
antiquark {QQ) potential. This is a typical pure gauge 
quantity (it is related to the expectation values of Wilson 
loops) which represents the starting point for many ap¬ 
proaches to the study of heavy quarkonia, tipically within 
a non-relativistic approximation. The static potential 
is usually expressed in terms of the so called Cornell 
parametrization [l^ : 

V(r) = ar . (1) 

r 

where a is the Coulomb coupling and a the string tension. 

Such a non-relativistic approach is reliable only for 
heavy quark bound states, for which the interaction en¬ 
ergy is not a large fraction of the total mass, and for 
not too large magnetic fields. As a rule of thumb one 
can use eBh/{m?c^) <C 1 (where m is the heavy quark 
mass) to estimate which magnetic fields can be explored 
in this approach. In our worst case (charm quark and 
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\e\B = 0.3 GeV^) the previous ratio is about 0.18, 
a numerical value which is roughly equal to the ratio 
{Mj/q, — 2mc)/Mj/q,, which gives an estimate of the va¬ 
lidity of the non-relativistic approach also for i? = 0. For 
bottomonia states eBh /gets smaller by more than 
one order of magnitude. 


II. THE QUANTUM TWO-BODY PROBLEM IN 
EXTERNAL MAGNETIC FIELD 

The non-relativistic Hamiltonian of two particles of 
masses rrii and charges qi (f = 1, 2) in an external mag¬ 
netic field is: 


The magnetic background field breaks rotational in¬ 
variance, hence one may expect the emergence of 
anisotropies in the potential, which is central otherwise. 
The issue is clearly related to quark loop effects, since 
gluons are not directly coupled to the magnetic field, and 
has been investigated within various model studies [Ti¬ 
ll^: the Coulomb coupling has been predicted to change 
in the transverse direction llR 12211 . in the longitudinal 
direction or in both |25l|: regarding the string 

tension, while string theory studies do not predict an 
influence of the magnetic field on it [26l - l^ . other ap¬ 
proaches do (see, e.g.. Refs. [2l|, |2^). 

Further insight into the question has been provided 
by a recent lattice QCD investigation [3l|: the potential 
gets steeper in the directions transverse to the magnetic 
field and flatter in the longitudinal one. In particular one 
observes a larger (smaller) string tension in the transverse 
(longitudinal) direction, while the opposite happens for 
the Coulomb coupling [^ . 

The purpose of this study is to reconsider the computa¬ 
tion of the heavy quarkonia spectrum and mixings in the 
presence of a magnetic field, in the light of the existing 
anisotropy in the static quark-antiquark potential. Since 
standard magnetic field effects decrease with the mass of 
the valence quarks, while the static potential remains un¬ 
changed, we might expect that corrections be relatively 
more important for bottomonia than for charmonia. 

We will follow the standard non-relativistic two-body 
approach (see, e.g.. Refs. [H, for reviews on this 
subject), which can be tackled using standard numeri¬ 
cal methods. In par ticular, our treatment will be close 
to that of Ref. [l^, apart from the presence of the 
anisotropy, which will be included according to the lat¬ 
tice results of Ref. [3l|. Since in Ref. [3l| the potential 
has been studied only along two directions (transverse 
and longitudinal), we shall adopt the simplest possible 
parametrization which takes into account generic angles. 


H = , 2 ) 

2 = 1 ^ ^ 

+ V {xi,X2) - (Ml + H2)- B , 

where Mij magnetic moments of the two 

particles. The presence of a vector potential A{x) 
makes the system not invariant under translations, i.e. 
Xi —>■ Xi + OL. Even in the case of an uniform mag¬ 
netic field, both the canonical and the kinetic momen¬ 
tum P = ~ ) do not commute with the 

Hamiltonian in Eq. ©■ 

An invariance group of H is obtained by simultane¬ 
ously performing a coordinate translation and a gauge 
transformation, i.e. by changing both Xi and The 
generator of this transformation is the pseudomomentum 
operator (see Q for more details) and can be written 
in a particularly simple form in the symmetric gauge^ 
A{x) = X a:, where it is given by: 

K X x^^ . ( 3 ) 

From this expression, it is not obvious that different com¬ 
ponents of the pseudomomentum are commuting observ¬ 
ables, and in fact this is not in general true. This is 
however what happens for globally neutral system (in 
particular for the QQ system to be studied in this pa¬ 
per), since the following commutation relations hold 

[Kj,ki] = -i{qi + q2)f-jtkBk ■ ( 4 ) 

From now on we will explicitly restrict to the case of a 
particle-antiparticle system and we will adopt the nota¬ 
tion q = qi = —q2 for the electric charge and m for the 
particle masses (see e.g. Ref. [l^ for the general case). 
Written in terms of the relative coordinate r = Xi — X 2 , 
the pseudomomentum takes the form 

k = Pi + P2 + ^qB xr = P + qBxr, ( 5 ) 


The paper is organized as follows. In Sec.|H]we will re¬ 
view the non-relativistic approach to the two-body prob¬ 
lem in an external magnetic field. In Sec. IHII we will 
discuss the form of the static potential in a magnetic 
background and our parametrization for its angular de¬ 
pendence. In Sec. m we will discuss our numerical ap¬ 
proach to the problem (more details and results from a 
test on the harmonic oscillator are reported in the Ap¬ 
pendix) and present results for charmonia and bottomo¬ 
nia. Finally, in Sec. |Vl we will draw our conclusions. 


From this expression it appears natural to adopt the fol¬ 
lowing ansatz for the eigenstate $ of the two body Hamil¬ 
tonian: 


<i>(i?, r, cr) = exp 


K — 2^B X r 


R 


4'(r, cr) 


( 6 ) 


^ See [ 35 II for the form of the pseudomomentum in a generic gauge. 
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where K denotes the eigenvalue of the operator R 
is the coordinate of the center of mass and cr is a short¬ 
hand for the spin variables. Using the expression Eq. m, 
the Hamiltonian in Eq. m can be rewritten in the form 
(acting on the reduced wave function 'I'): 


H =— - -^{K X B) ■ r - — + 

2M jvr ’ 2fi 

+ X r)2 + Vir)- (/x^ + ■ B , 


( 7 ) 


where M = 2m and n = m/2 are respectively the total 
and the reduced mass of the two body system. Because of 
the lack of translation invariance, in this Hamiltonian the 
center of mass and the relative motion are not decoupled, 
unless K X B = 0. 

Given this nontrivial dependence on the K value, one 
needs a prescription to subtract the energy associated 
with the center of mass motion. We will follow the idea 
presented in Ref. [l^ : given the eigenvalue E of Eq. ©, 

/ f >\2 

we will consider E — as the corrected binding energy, 
with the kinetic momentum expectation value being com¬ 
puted by using Eq. ©• 

Einally, let us discuss the role of the spin interaction 
term in Eq. ©• Quark magnetic moments may be ex¬ 
pressed as /Xj = where g = 2 is the quark g-factor 

in the non-relativistic approximation, /x^ = qi/2mi is the 
quark magneton and Si = cr®/2 is its spin. Therefore, in 
the case of a quark-antiquark pair 


- {fii+fi2) ■ B =--^{cT^ - cr^) ■ B . ( 8 ) 

It can be verified that this term induces a mixing between 
the singlet |00) and the triplet |10) spin states. Indeed, 
in the presence of a magnetic field B — Bz, the operator 
in Eq. (|5]) has the only non-zero matrix element between 
spin states given by 

(00|(/Xi+/x,).B|10) = -^ . (9) 

Such a term increases the energy of the triplet state and 
decreases that of the singlet. When Eq. m is used for 
a noncentral potential, or when a spin-spin interaction 
is also present, some caution is needed, since possible 
mixings between different orbital states can be present 
(see Appendix. for more details). 


III. THE QQ POTENTIAL 

It is well known that the main features of the quarko- 
nium spectrum can be understood by using a central po¬ 
tential between the heavy quarks of the so called Cornell 
form reported in Eq. ©. This parametrization was in¬ 
troduced in Ref. [l^ on phenomenological bases, as the 
simplest form of the potential that takes into account 
both the short distance perturbative contribution and 


QQ 

State 

Name 

Mass [MeV] 

cc 

us'd 

be 

2980.3 ± 1.2 


CSi 

J/P 

3096.916 ±0.011 


l^Po 

XcO 

3414.75 ±0.31 


l®Pi 

Xcl 

3510.66 ±0.07 


1^P2 

Xc2 

3556.20 ± 0.09 


l^Pi 

he 

3525.38 ±0.11 


2^ So 

gc{2S) 

3639.4 ± 1.3 


2^ Si 

P{2S) 

3686.109 ± 0.02 

bb 

l^So 

gb 

9390.9 ± 2.8 


CSi 

T 

9460.30 ± 0.26 


l^Po 

XbO 

9859.44 ± 0.74 


l^Pi 

Xbl 

9892.78 ± 0.57 


1^P2 

Xb2 

9912.21 ±0.57 


l^Pi 

hi 

9898.3 ± 1.1 


2^ So 

76(25') 

9999.0 ± 3.5 


2^ Si 

T(25) 

10232.26 ± 0.31 


TABLE I: Lowest part of charmonium and bottomonium mass 
spectrum. Data from [3^ . 


color confinement. It was later realized that the form 
Eq. ([T]) of the potential correctly describes the spin aver¬ 
aged potential between two static (i.e. infinitely massive) 
quarks as estimated from the first principles of QCD by 
its lattice formulation (see, e.g.. Ref. 

To properly describe the fine structure of the quarko- 
nium spectrum further spin dependent terms have to be 
added to the potential. All such corrections to the static 
quark potential can in principle be obtaine d by an expan¬ 
sion in the inverse quark mass (see Ref. |33[), however 
their precise functional form is generally not well known 
and difficult to extract from LQCD computations. Be¬ 
cause of this lacking of precise theoretical information, 
several different parametrization of these terms exist, 
whose coefficients are fixed by comparing with experi¬ 
mental results. 

Since we are mainly interested in the lowest quarko- 
nium levels (see Tab.lJ), the most important contribution 
to be added to the Cornell potential is the spin-spin in¬ 
teraction, responsible e.g., for the mass splitting between 
the J/'i> and rjc levels of charmonium. We will adopt the 
parametrization 


Uto- = (cTl • 0-2)76 


( 10 ) 


where P and 7 are phenomenological constants, whose 
values are reported in Tab. nil This form of the spin-spin 
interaction is supported by the lattice result presented in 
Ref. [131 and was previously used in Ref. [T^ to study 
the influence of a magnetic field on the meson spectrum. 

While in previous studies the static Cornell potential 
in Eq. © was assumed to be independent of the magnetic 
field, in this work we will consider also the anisotropy ob¬ 
served in Ref. 311. However, since in Ref. [llj only the 
values of the potential along the coordinate axes have 
been investigated, we are forced to make an ansatz on 
the form of the potential, with the constraint that it re¬ 
produces the observed behaviour on the axes. To this 










4 


0 

J^Oxy 

^Oxy 

AOz 

COx 

a 

a 

0.29T0.02 

-0.24T0.04 

0 0 

0 0 

-0.34T0.01 

0.24±0.03 

1.5T0.1 

1.7±0.4 


TABLE II: Coefficients A, C estimated in [3l| . 


aim, we propose the following anisotropic form of the 
static potential terms {B is directed along the z axis) 


a a 



where the scaling parameters exy{B) and ez{B) are func¬ 
tions of the external magnetic field. Such a parametriza- 
tion is inspired by that of the electrostatic interaction 
in the presence of an anisotropic dielectric constant (see 
Ref. @ §13). It is possible to rewrite the expression 
above in a more usual fashion, by absorbing the angular 
and B dependences into the a and a parameters 


V(r,9,B) = -^^^ + aie,B)T , ( 12 ) 


where 9 is the spherical azimuthal angle and 
a{9,B) = 


>yi + e(“Hi?)sin2(0) 
(7{9,B) = crej'"^ {B)\Jl + {B) sin^ (9) 


(13) 


with related to the scaling parameters previously de¬ 
fined by setting ei = y/e^ and €2 = e^y/cz — 1. According 
to Ref. , when increasing the magnetic field the string 
tension gets larger in the transverse directions, while it 
decreases in the longitudinal z direction; the opposite be¬ 
haviour is observed for the coefficient a of the Coulomb 
term. To reproduce such a behaviour, the coefficients 
must be positive. In particular, they are related to 
the coefficients A, C estimated in Ref. ^ (see Tab. [II]) 
by the relations 

-1 


{^\b) = 

(l + A“«(|e|P)^“') 

PiB) = 

/ I + A“-(|e|P)'^“" 


[^\b) = 

(l + A'^^(|e|P)^’^) 

■^2\B) = 


1, I + A-K|e|P)C- 


- 1 


- 1 . 


(14) 


A graphical representation of the anisotropic potential is 
shown in Fig. [I] 

Finally, let us briefly discuss the role of the spin-spin 
interaction term (HQl) within the anisotropic potential pic¬ 
ture. In principle, both parameters, /3 and 7 , may ac¬ 
quire a dependence on the external magnetic field like 



FIG. 1: Contour map of the potential in Eq. d on the X — z 
plane for B = 0.6 GeV^, with B directed along the z axis. 
Parameters used are the ones for charmonium (see Tab. EH). 


that observed for the string tension and the Coulomb 
coupling. Such a dependence has not yet been investi¬ 
gated and no ansatz is available. However, since this 
term represents a spin-dependent relativistic correction 
of order [H, Hj, IH, li^, it is reasonable to expect 
that also the B-dependence be of the same order in the 
quark mass m, hence much weaker than that induced on 
the spin-independent part of the potential. Therefore in 
the following we will make use of the same values for 
P and 7 adopted in Ref. [T^- A similar attitude will be 
adopted with respect to other relativistic corrections, like 
the spin-orbit term. 


IV. COMPUTATION AND NUMERICAL 
RESULTS 

In this section we report our results for the dependence 
on the magnetic field of the masses of the IB and IP 
states for both charmonium and bottomonium. 

To this purpose we used the Hamiltonian in Eq. 0, 
with a magnetic field B = Bz and pseudo-momentum 
K = Kx, with the anisotropic potential parametrized 
as in Eq. m, together with the spin-spin interaction 
term Eq. (uni). Values of the parameters adopted are 
reported in Tab. ED and have been fixed according to 
Ref. [ 13 , thus enabling a direct comparison with the case 
in which the anisotropy of the potential is neglected. Nu¬ 
merical values of the anisotropy parameters, defined in 
Eq. (fHl) . have been fixed by using the coefficients shown 
in Tab. El To take into account the uncertainty associ- 
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13 

a 

a 

m 

cc 

2.060 GeV 

1.982 GeV 

0.312 

0.174 GeV^ 

1.29 GeV 

hb 

0.318 GeV 

1.982 GeV 

0.421 

0.210 GeV^ 

4.70 GeV 


TABLE III: Parameters used for the potential (same as in 
Ref. HI]). 

ated with these parameters, several simulations were per¬ 
formed, corresponding to different combinations of these 
values. 

In order to allow for a better physical comprehension 
of the results, the mass spectrum will be shown as a func¬ 
tion of the mean kinetic momentum (P), instead of us¬ 
ing the pseudo-momentum K. To this purpose we solved 
the system for several values of iT, computed the value 
of the mean kinetic momentum for each eigenstate and 
interpolated these values to the desired point. Details 
of the numerical algorithm used to extract eigenstates 
are reported in Appendix 1X1 together with a test of the 
algorithm in an analytically solvable case. 


A. Charmonium states 

We start by studying the r]c and J/ip states, which cor¬ 
respond to the IS' states of the charmonium at vanishing 
magnetic field. The experimental masses of these states 
are well reproduced by using the potential parameter re¬ 
ported in Tab. IIIII and their behaviour as a function of 
the external magnetic field B is shown in Fig. [5] The spin 
components Sz = ±1 of the J/ip stay degenerate also for 
B ^ 0 but they split from the = 0 component, which 
mixes with the state. 

In order to isolate the contribution of the anisotropy in 
the potential, we reported in Fig.[2]also the masses com¬ 
puted without taking into account the anisotropy (i.e. 
with = 1 and £ 3 “^ = €^ 2 '^ =0). It can be 

seen that in all cases the masses are increased by taking 
into account the anisotropy; as a consequence, for the qc 
the dependence of the mass on B gets milder, while the 
opposite effect is observed for the J/ijj states. We note 
that for the qc and the Sz = ±1 components of the J/ip 
the effect induced by the anisotropy is of the same order 
of magnitude as the effect expected when no anisotropy 
is present. 

A non-vanishing kinetic momentum has an effect anal¬ 
ogous to that of the anisotropy: it reduces the depen¬ 
dence on B of the qc mass while it increases that of the 
J/'0 states. This effect is however quite small for the val¬ 
ues of the kinetic momenta explored and it is almost of 
the same magnitude of the uncertainties associated with 
the anisotropies. It cannot be excluded that this effect 
increases for larger momenta, but to systematically ex¬ 
plore this regime a fully relativistic treatment would be 
required. 

As previously noted, the magnetic field introduces a 
mixing between the qc state and the Sz = 0 component 


(a) 



(b) 



(c) 



FIG. 2: Behaviour of the masses of the charmonium IS states 
with respect to the magnetic field, both with and without the 
magnetic anisotropy in the potential and for two different val¬ 
ues of {P). Data points represent the mass extracted making 
use of the central values of the parameters in Tab. m while 
shaded regions take into account uncertainties on the param¬ 
eters. (a) The qc singlet state, (b) The Sz = ± 1 components 
of the J/i/> state, (c) The Sz = 0 component of the J/'tp state. 


of the J/ip state. The behaviour of this mixing with re¬ 
spect to the magnetic field and the pseudomomentum K 
is shown in Fig.|31 while the dependence on the magnetic 
field is quite strong, the mixing turn out to be almost in¬ 
dependent of the value of the pseudomentum. Moreover, 
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contrary to what happens for the masses, the mixing is 
almost insensitive to the anisotropy of the potential. 



|e|B [GeV2] 

FIG. 3: Square modulus of the mixing between the Sz = 
0 component of the J/i/’ state and the rjc state, both with 
and without the magnetic anisotropy in the potential and for 
several values of the pseudomomentum K. 



FIG. 4: Behaviour of the masses of the charmonium IP ex¬ 
cited states with respect to the magnetic field. We report for 
simplicity only data obtained by using the anisotropic form of 
the potential, the general features of this picture are present 
also in the standard case of the central potential. 

We now present some preliminary data for the IP 
states. In this case the fine spectrum of known IP char¬ 
monium levels is not correctly reproduced for P = 0, a 
fact that is likely due to the absence of the spin-orbit 
coupling and other relativistic corrections in the used 
Hamilonian. If we assume that these relativistic terms 
do not depend on the external magnetic field, the com¬ 
puted variation of the mass Ato(P) due to a nonvanish¬ 
ing B can be used to shift the known B = 0 values of 
the masses to obtain the spectrum at P ^ 0. It is clear 
that this is an approximate procedure and that a more 
precise study will be required to quantify the systematic 
error introduced in this way; nevertheless this method is 
expected to give reliable informations in the limit of large 
quark mass m, since relativistic effects are suppressed by 
inverse powers of m. The results obtained in this way are 


shown in Fig. |4] and the main difference with respect to 
the 15 case is that in the IP case several level crossings 
happen by increasing the magnetic field. From Figs. 0] 
and 0] it can be seen that the gap between the h states 
and the J/tp states is strongly reduced by increasing the 
magnetic field, a fact that may have significant conse¬ 
quences on branching ratios of the h meson. 

B. Bottomonium states 

The effect of the anisotropy in the potential on bot¬ 
tomonium states is qualitatively similar to that on char¬ 
monium states, however its quantitative relevance is 
much larger, as can be seen from the spectrum of the 15 
states in Fig. [SJ the mass shift due to the anisotropy is 
about one order of magnitude larger than the one due to 
quark magnetic momenta only. This is not unexpected, 
since the effect of the anisotropy in the potential is inde¬ 
pendent of the quark masses, while quark magnetic mo¬ 
menta go to zero like 1/m. As a consequence the relative 
effect of the anisotropy is much stronger for bottomonium 
than for charmonium. 

In Fig.[6]we show the behaviour of the mixing between 
the singlet state ryf, and the = 0 component of the T 
as a function of the magnetic field. Although the value 
of the mixing is reduced by about a factor 3 (~ mb/me) 
with respect to the charmonium case, its value is still 
largely independent on the anisotropy. 

A preliminary analysis of the IP states is shown in 
Fig. [71 it is obtained by the same procedure adopted for 
the charmonium, which for bottomonium is expected to 
be more reliable. Also for bottomonium some level cross¬ 
ings can be seen, however in this case the uncertainties, 
due to the propagation of the error on the anisotropy, 
become quite large with respect to the typical level sep¬ 
arations. 


V. SUMMARY AND CONCLUSIONS 

The main purpose of this work has been the study of 
the possible influence of strong external magnetic fields 
on the heavy-meson mass spectrum. Following the ap¬ 
proach of Ref. [l^l, we used a non-relativistic model to 
extract the lowest lying cc and bb states and the magnetic 
field induced mixings among them. As a new ingredi¬ 
ent, we considered the anisotropies induced by the mag¬ 
netic field at a non-perturbative level in the static quark- 
antiquark potential, as determined by a recent LQCD 
investigation [3l|. Moreover, we extended our analysis 
to the first excited IP states. 

Various approximations are involved in our computa¬ 
tion: relativistic corrections were only partially included 
and any possible dependence of their magnitude on B 
was neglected. Moreover, the exact angular dependence 
of the static potential was obtained by a suitable inter¬ 
polation of data reported in Ref. [3l|. 





7 


(a) 



CM 


O 


c* 

V 



(b) 



|e|B [GeV^] 

(c) 



FIG. 5: Behaviour of the masses of the bottomonium IS 
states with respect to the magnetic field, both with and with¬ 
out the magnetic anisotropy in the potential and for several 
values of (P). (a) The rjb singlet state, (b) The Sz = ±1 
components of the T state, (c) The Sz = 0 component T 
state. 


A first issue of phenomenological interest regards the 
mixing between different states, since that should lead 
to a modification of both decay patterns and production 
rates. For instance, results suggest a contamination of 
the semi-leptonic decay channels between the r]c and the 
J/4> and between the ryt, and T mesons. However, in this 
case such mixings are observed even without taking into 


FIG. 6: Square modulus of the mixing between the Sz = 0 
component of the T state and the T]h state, both with and 
without the magnetic anisotropy in the potential and for sev¬ 
eral values of the pseudomomentum K. 



FIG. 7: Behaviour of the masses of the bottomoniumm IP 
excited states with respect to the magnetic field. 

account the potential anisotropy m , and we have veri¬ 
fied that its inclusion does not lead to significant quanti¬ 
tative changes. 

On the other hand, results suggest that the presence 
of the anisotropy in the potential has significant effects 
on the heavy-meson mass spectrum. An increase of the 
masses is generally observed, with respect to the case in 
which the anisotropy is not included, the increase being 
of the order of 30 — 40 MeV at the maximum magnetic 
field explored, \e\B ~ 0.3 GeV^. The effect is more dra¬ 
matic for bottomonium states, where it can even change, 
in some cases, the sign of the dependence of the mass 
on B: the reason is that standard magnetic interaction 
terms in the Hamiltonian of the system (see Eq. ©) are 
suppressed by the inverse of the quark mass, while the 
non-perturbative effect on the static potential is mass in¬ 
dependent, hence it becomes dominant for bottomonium. 

Finally, a new effect pointed out by our study regards 
the possible crossings of IP states as a function of B. 
Such crossings are observed even without taking into ac¬ 
count the anisotropy of the potential, however its pres¬ 
ence makes them clearer. We stress that our results for 
















IP states are still preliminary: in particular, regarding 
the level crossings, we assumed the experimentally ob¬ 
served spectrum at P = 0, which can be reproduced for 
such states only by taking into account relativistic cor¬ 
rections like the spin-orbit one. On the other hand, the 
spin-orbit term was not taken into account to compute 
the P-dependence of the spectrum, therefore a significant 
P-induced correction to the this term could change the 
scenario. 

One should consider that present results are obtained 
using the T = 0 form of the potential and assuming a con¬ 
stant and uniform magnetic field. The use of the T — 0 
potential is justified only for hard processes taking place 
in the initial stages of non-central collisions, before the 
strongly interacting medium thermalizes. Future lattice 
simulations could provide information on magnetic field 
effects at finite T, which could be relevant to produc¬ 
tion and decay rates of heavy quark states in the ther- 
malized medium; instead, it will be relatively more diffi¬ 
cult to take into account the possible effects of inhomo¬ 
geneities or time dependence of the magnetic field dis¬ 
tribution. Further improvement on present results could 
also be obtained once lattice simulations provide infor¬ 
mation about the exact angular dependence of the static 
potential (which in the present study was partially based 
on an ansatz) and, possibly, about the spin-dependent 
part of the potential. 


(and the corresponding among those correspond¬ 
ing to states such that ($a|C) 0- By using an initial 

wave function C(r) with specific symmetries we can thus 
select the state we are interested in and, in this work, 
hydrogen-like wave functions were used to this purpose 
when B — 0. States for B ^ 0 were extracted by adiabat- 
ically switching on the magnetic field in an initial stage 
of the evolution. This procedure turned out to work well 
for the low lying states, however numerical instabilities 
emerge when applying this technique to higher excited 
states. 

From the numerical point of view the evolution is per¬ 
formed by introducing a temporal and a spatial lattice 
spacing (denoted by dr and a respectively) and by ap¬ 
proximating derivatives in Eq. jAT]) with finite differ¬ 
ences. In this way variables at time r + dr can be written 
in term of the ones at time r (see [4l| for details). 

After every time evolution step we get an estimate of 
the bound state energy and of other properties of the 
state we want to study; the process ends when the vari¬ 
ation of these quantities during a time interval Ts goes 
below a fixed threshold. The time interval Tg was set to 
1/M, where M is a rough estimate of the mass for the 
state we are interested in. The precision of 1 MeV was 
used as a stopping criterion for the energy determination 
(this value can be used as the uncertainty in the energy 
computation). 
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Appendix A: The numerical algorithm 

The algorithm used to numerically compute eigenval¬ 
ues and eigenfunctions of the reduced Hamiltonian © is 
the Finite Difference Time Domain method (FDTD) de¬ 
scribed, e.g., in [Uld^. The main idea of this approach 
is the following: once the Schrodinger equation is Wick 
rotated to imaginary time 

(J-+H^^ij{r,T)=0 , (Al) 


9.36 

> 

U 9-28 
w 9.24 
9.2 



0 0.1 0.2 0.3 0.4 0.5 

a [GeV'l] 


9.34 
^ 9.32 

I 9-3 

w 9.28 
9.26 



0 0.002 0.004 0.006 0.008 

dt [GeV'l] 


9.34 

? 9-32 

S 9.3 
“ 9.28 
9.26 



0 0.1 0.2 0.3 0.4 0.5 

a [GeV-l] 


the formal solution of the problem with initial condition 
'I'(r, 0) = ({r) is 

vl-(r,r) = ^($„|C)$„(r')e-^“" , (A2) 

Ct 

where and are the eigenvalues and eigenfunctions 
of the Hamiltonian. By looking at the large time be¬ 
haviour of 4'(r,T) we can thus identify the lowest Cq, 


FIG. 8: Test of the continuum limit for parameters uj = 1.5 
GeV, |e|i3 = 0.9 GeV^ and K = 5.0 GeV. Data correspond to 
the energy of the eigenstate , while dashed lines are fit 
curves, (a) First dr —^ 0 and then a —>■ 0, the final continuum 
value is Ea = 9.3205(13) GeV. (b) First a —>■ 0 and then 
dr 0, El, — 9.3206(13) GeV. (c) a — >■ 0 with a oc dr^, 
Ec = 9.3201(12) GeV. 
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(a) 



£ 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

|e|B [GeV^] 


bital momenta, however, turns out to be rather small: 
in all the studied cases the overlap = 0|£ = 1) was 
compatible with zero within machine precision and the 
overlap (£ = 0|^ = 2) was at most of the percent level. 

All the results presented in the main text where 
obtained by using a lattice of physical spatial extent 
V = (30 GeV“^)^ ~ (6 fm)^ and several values of 
the lattice spacings, in order to extract the contin¬ 
uum limit. The spatial lattice spacings used were a = 
0.250,0.375,0.500,0.625 GeV“^ and the temporal lattice 
spacing was fixed by the relation dr = ma^/20, m being 
the quark mass. 


(b) 


1. Test of the algorithm 



4_6 u_^^^^^^^^^^_I_^^^^^^^_ 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

|e|B [GeV^] 


FIG. 9: Behaviour of the masses for the eigenstates $( 0 , 3 ) with 
respect to the magnetic field. Dots correspond to numerical 
data extrapolated to the continuum, dashed curves represent 
the analytic expected behaviour, (a) Data corresponding to 
\K\ = 0. (b) Data corresponding to \K\ = 5 GeV, K = IJfl®. 


The numerical algorithm just described concerns the 
orbital part of the Hamiltonian and can be applied, once 
the total spin is fixed, to the reduced Hamiltonian Eq. © 
(without the ■ B terms) together with the V^a potential 
defined by Eq. cni). By using these eigenstates we can 
then evaluate the matrix elements of the operator (/x^ -|- 
( 12 ) B and diagonalize this matrix to extract the final 
eigenvalues and eigenstates. 

The operator (/x^ + M 2 ) ' non-vanishing matrix 

elements only between the spin states |00) and |10) (see 
Ref. 0 for details), however some care is needed here 
since, due to the specific form of adopted, noth¬ 
ing prevents a mixing between states with different or¬ 
bital momenta to occur (this problem is relevant also for 
B = 0). The effect of this mixing between different or- 


We report here the results of some tests performed to 
verify the correctness of the algorithm implementation. 
As a testbed, we used the case of the harmonic potential 

V{r) = , (A3) 

which is analytically solvable (see Ref. 0111) and thus 
allows for a direct check of the numerical data. Following 
the notation of Ref. 0 ] , eigenstates are classified by the 
quantum numbers K, n±, and Sj,, where K is the 

pseudo-momentum defined in Eq. s is the total spin 
of the system, its projection on the z-axis and n_L, n^, I 
are quantum numbers specific for the harmonic potential. 
For fixed spin s, the ground state will be denoted by 
and corresponds to the case n_L = = £ = 0 

(the superscripts ± and 0 will be used to denote the Sz 
components). 

As a first test we verified that, for fixed physical pa¬ 
rameters, the same result is obtained by performing the 
continuum limit in different ways: it is possible to first 
extract the limit for a —>■ 0 at fixed dr and then send 
dr to zerol; alternatively, it is possible to exchange the 
order of the limits or to perform both limits together. 
Since in our implementation discretization errors are lin¬ 
ear in dr and quadratic in a, it is convenient to perform 
the limit by using a relation of the form dr oc c? with a 
fixed proportionality factor. 

The three possible ways to continuum extrapolate are 
compared in Fig. [8] and they nicely agree with each other 
within errors. In Fig. [S] the final results (i.e. continuum 
extrapolated and taking into account spin mixing) for 
the energies of the low lying states are compared with 
the theoretical expectations, and in all cases a perfect 
agreement is found. 
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